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Abstract 

We show that in 1+1 dimensional conformal field theories, exciting a state with a 
local operator increases the Renyi entanglement entropies by a constant which is the 
same for every member of the conformal family. Hence, it is an intrinsic parameter that 
characterises local operators from the perspective of quantum entanglement. In rational 
conformal field theories this constant corresponds to the logarithm of the quantum 
dimension of the primary operator. We provide several detailed examples for the second 
Renyi entropies and a general derivation. 
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1 Introduction 

Conformal field theories (CFT) in two dimensions play a very important role in the under¬ 
standing of quantum entanglement (see [1] for review). For example, in the vacuum state, 
one can compute entanglement measures like Renyi entropies analytically for an arbitrary 
interval. They are universal and determined by the central charge c and in particular the von- 
Neumann entropy has the famous | logL scaling with the size of the interval L. This relation 
has been tested quite extensively and serves as an efficient way to numerically determine the 
central charge for the CFT that governs particular critical points. 

The next natural step is to explore the entanglement in excited states. A particularly 
useful protocol to study that, is provided by global or local quenches [2j where one can 
follow the time evolution of entanglement measures in a state that differs from an eigenstate 
globally or locally respectively. Under such evolution, entanglement grows in the system and 
it is hard to study numerically, nevertheless, the power of (boundary) CFT allows to extract 
the universal features like the speed of the growth of entanglement. 
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Recently, inspired by the quench setup, another class of excited states by local operators 
was proposed in j3jj] They can be thought of as a milder version of a local quench and the 
growth of entanglement is less drastic. In fact, at late times, the Renyi entanglement entropies 
saturate to constants that can be used to characterize local operators by the way they change 
quantum entanglement in a given state. Detailed reviews can be found in pm and further 
results and applications to finite temperature and holography in [8J 0 [101, HU H2J 13], 

For primary operators in 1+1 dimensional rational CFTs, the growth of the entanglement 
of an interval can be computed analytically and it was proved to be the logarithm of the 
quantum dimension of the local operator PJ. Since in this setup we can uniquely decompose 
an operator into left and right movers (chiral, anti-chiral parts) this constant is equivalent the 
the left-right entanglement in this excited state. Intuitively, this is explained by the quasi¬ 
particle picture for the propagation of entanglement where the excitation can be thought of 
as insertion of the EPR like pair. Quasiparticles propagate in opposite directions and once 
one of them is in the entangling region it increases the entanglement with the rest by the 
entanglement of the pair. This picture is supported by the one point functions of the energy 
in locally excited state [H] . 

In this work we also focus on states excited by local operators in 1+1 d CFTs with 
Virasoro symmetry, and in the search for universal features of entanglement we ask how the 
above picture changes for descendant operators. As we will show, one-point functions of the 
energy behave differently in the states excited by descendants. Moreover, their conformal 
transformations, as well as their correlators, are much more complicated than for primaries 
and the appearance of the quantum dimension is not obvious. Nevertheless, we will show for 
descendants up to level 2 and provide a general argument that the constant contribution to 
the Renyi entanglement entropies is the same for the entire conformal family. 

This paper is organized as follows. In section 2 we briefly review the replica method with 
local operators in 2d CFT |l 9] and provide an example that illustrates the CFT technology. 
In section 3 we compute one-point functions of the energy density and the increase in the 
second Renyi entropy for descendant operators up to level 2. In section 4 we outline a general 
argument for why the contribution to the Renyi entanglement is a characteristic of the entire 
conformal family. In section 5 conclude and discuss the Schmidt decomposition and possible 
connection to topological entanglement entropy. Several details of the analysis are moved to 
the appendices. 

2 Replica trick for locally excited states 

The entanglement entropy of a subsystem quantifies the amount of information that we 
would forfeit if we were to loose access to the rest of the system. Imagine that the system 
of interest is in a pure state +) and suppose that we wish to quantify the entanglement 

1 see also 00 for a different setup with local operators 
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between a subsystem A and its complement A. As a first step we find the reduced density 
matrix obtained by tracing out the degrees of freedom in A 


Pa = Tr^J| . (2.1) 

Once we have constructed this matrix, we notice that if there is any entanglement between 
the degrees of freedom in A and those in A the system appears, to an observer having access 
only to A, to be in a mixed state. If that is the case, then the Von Neumann entropy of pa 

S A = -Tr (p A log p A ) , (2.2) 


is non-vanishing. We refer to this quantity as the entanglement entropy of A. In practice, 
we use the so-called the replica trick [I] and compute instead the Renyi entropies 

s?’= T W°g' Iv ra • (2.3) 


The entanglement entropy (2.2) can be extracted from the above expression by taking the 
n —> 1 limit. These Renyi entropies provide interesting measures of entanglement on their 
own right e.g. the min entropy S ^°°' ) and the purity [15]. 


It is natural to wonder what would be the effect on the entanglement between different 
parts of a system if we were to perturb it in some way. In this work we focus on the case 
where these perturbation is due to the insertion of a local operator. Recently, the replica 
method to compute the Renyi entropies has been generalized to deal with these scenarios. 
Below we summarize the relevant formulas and refer the reader for more details to mm 


Consider the ground state of a 1+1 dimensional CFT on a line. We split the space into a 
finite interval A = (7 1; l 2 \ of length L = l 2 — l\ and its complement. Then at t = 0 we insert 
a local operator O at x = — l and let the system evolve. Without losing generality, we can 
shift the interval to A = [0, L\. The resulting density matrix is given by 


p(t,l,L,e ) = M • e- iHt e- eH O{0, -l) |0) (0|O + (0, -l)e~ tH e iHt 

= Af ■ 0(w 2 ,w 2 )\0)(0\0^(wi,w l ) . (2.4) 

where e 1 plays the role of an UV regulatoi0 and J\T is a normalization that ensures that 
Tr p — 1 . In the second line the insertion points are defined as 


wi = i(e — it) — l, w 2 = —i(e + it) — l, 

Wi = —i(e — it) — l, w 2 = i(e + it) — l. (2.5) 

Hereafter, we omit the dependence on /, L, e and write simply p(t). 

Analogously to the local quench setup, we assume that all length and time scales are much larger than 
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The increase in the n-tli Renyi entanglement entropy of the interval A wrt to the ground 
state due to the local operator is given by [3] 


A S ( A n) = 


n 


log 


(Q(w 1 ,Wi)Oi(w2,W 2 ) ■ ■ ■ 0(w2n-liW2n-l)0 ] (w 2 niW2n))l 

({^(Wi, wi)0(w 2 , W 2 ))y, 1 ) n 


( 2 . 6 ) 


where E n is the n-sheeted surface with cuts on each copy corresponding to A, and Eh is a 
single copy with an interval cut A. Thus, we are faced with the task of evaluating a 2n-point 
function on E n . As shown in [9], this can be computed by using the uniformization map 

= ~~y , ( 2 - 7 ) 

that takes £ n to the complex plane. Notice that the transformation properties of O under 
the above map are expected to play an important role. 

Let us focus on the change in the second Renyi entropy A S K A ; . In this case we are expected 
to evaluate the operator’s four-point function on a two-sheeted surface. We use (2.7) with 
n = 2 to map E 2 to the plane, where conformal symmetry ensures that these correlator 
can be written in terms of the cross-ratios of the four insertion points. Notice that under 
the uniformization map the insertion points (2.5) become Z 4 = —z 2 and z 3 = —Z\. From 
equations (2.5) and (2.7) and in the limit e —>• 0 the cross ratios are given by [9j 

L 2 e 2 L 2 e 2 


£12-34 

z =-~ 1 — 


z ~ 


-13-24 4 (l - t) 2 (L + l-t ) 2 ’ 4 (l + t) 2 {L + l + t ) 2 ’ 

provided that i £ [l, L + /]. On the other hand, if t (f [/, L + l\, then the cross-ratios read 

L 2 e 2 _ z 12 z 3 4 L 2 e 2 


( 2 . 8 ) 


z ~ 


4 (l -t) 2 (L + l -t) 2 ' 


z = 


-13-24 4(1 + t) 2 (L + l + t ) 2 


(2.9) 


Hence, there are two possibilities, (z,z) — (1,0) for t e [l,L + l\, while (z,z) — (0,0) 
otherwise. 

The authors of [9] showed that if O is a primary operator, then 

AS%\z,z) = -\og[\z(l-z)\ lh g(z,z)] . (2.10) 

where Q (z, z) is the canonical four point function of the operators on the plane. As discussed 
above, we must evaluate this result for the regimes (z,z) — (0,0) and (z,z) — (1,0). In the 
hrst regime there is no increase in entanglement 


asJ, 2) (o,o) = 0. 


( 2 . 11 ) 


Now, since the two regimes are related to one another by the modular transformation z —)• 
1 — z, in rational CFT, it follows that 


ASf(l,0) = - log F 0 „[O], 


( 2 . 12 ) 
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where F 0 q[O] is the component of fusion matrix for the identity [16]. Moreover, F 00 [O] 
corresponds to the inverse of the quantum dimension do of the operator H3. Therefore, one 
gets 


A c(2) _ f 0 if t £ [l, L + 1] 

' \ log do if t e [l,L + 1} 


(2.13) 


This behavior of the Renyi entropies has been interpreted in terms of a quasi-particle pic¬ 
ture where the insertion of the local operator corresponds to creation of the EPR-like pair 
whose members propagate in the opposite directions (to the left and right). A non-trivial 
contribution to the Renyi entropy comes from times when one of the members is inside the 
interval A while the other one is outside. 


In this work we will test these results and quasi-particle picture further by considering 
excitations by descendant operators. 


2.1 Example: EPR-primary 

In this section we consider a simple illustrative example of a primary operator that is a 
generalization of the operators studied in [H], which will also serve to check our results for 
descendants. Let us start with a state locally excited by the primary operator 

O a = Vd e iV * a * ± e-iyfa* ; (2.14) 


where 0 is a free massless scalar and a 
dimension of the operator is h = h = a 2 
EPR-primary. 


G [0,1]. It is easy to verify that the conformal 
. Hereafter, we refer to the operator (2.14) as the 


Replica method 

For the sake of simplicity we consider the second Renyi entropy. As discussed in the previous 

( 2 ) 

section, since O a is a primary operator the increase of S \ ; is given by 

AS^ 2) = - log 0(1 - z)\ 4h G a (z,z)] , (2.15) 

where Ga(z,z) is defined by the four-point correlator of the EPR operators on the complex 
plane via 

(O a (zi, z^Oliz*, z 2 )O a (z 3 , z 3 )O f a (z 4 , z 4 )) = \z 13 z 24 \~ 4h G a {z, z ), (2.16) 

or more explicitly we define 

G a (z,z) = (O a \O a (z,z)O a (l,l)\O a ) = \z(l-z)r 4h Fa(z,z). (2.17) 

The function T a (z,z), which encodes the non-singular part of G a (z,z) can be computed 
explicitly for the EPR operators (see Appendix [A]) and it is given by 

T a {z, z) = a 2 + (1 — a) 2 + 2a(l — a) (|^| 8h + |1 — z\ >h ) . (2.18) 
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(2.19) 


In terms of this function, the change in the Renyi entropy takes the compact form 

exp(-AS^ 2) ) = F a {z,z). 


Finally, for times t ^ [l,L +1\ we have both (z, z) — (0,0) and the change in the second 
Renyi entropy vanishes. On the other hand for t G [/, L + l\ the cross-ratios are approximately 
(z, z) ~ (1, 0) , which gives a non-trivial contribution to the second Renyi entropy 



0 

- log [ a 2 + (1 - a) 2 } 


t [l, L + l] 
t £ [/, L T Z] 


( 2 . 20 ) 


Notice that the non-vanishing contribution to the second Renyi is maximized by an EPR 
primary with a — 


Left-right entanglement 


In the following, we show how the same result can be obtained from the entanglement between 
left and right moving sectors in the state 


O a |0) L 0 10)* = 


Jy/2a<f) L \ 
/L 


D i\/2ci<t>R ^ 


R 


±VT^ 


0 -iy/2a(t> L \ 

/l 


D -iy/2a<t> R ^ 


R 


( 2 . 21 ) 


where we used the decomposition of the scalar field 4>(t,x) = 0x(t + x) + — x). Tracing 

out the right-moving sector we find that the reduced density matrix of the left-movers, which 
is simply 

Pl = diagja, 1 — a} . (2.22) 


Hence, the n-th Renyi entropy is given by 


S™ = 


1 


1 — n 


logTr p n L = 


1 — n 


log [a n + (1 - a)"], 


(2.23) 


which for n = 2 reproduces Eq. (2.20). 


Clearly the left-right decomposition of the excited state gives a quick way to find the 
reduced density matrix. Once we have itrl we can study the Renyi entropies for different 


values of n and see how they are affected by the local excitation (2.14). For instance the 
Hartley entropy (n = 0) is given by S[ Q> = log(2) and sets an a-independent upper bound 
for all the Renyi entropies. It corresponds to the log of the dimension of the reduced density 
matrix pl■ In turn, the entanglement entropy (n = 1) equals the binary entropy 


= H- 2 (a) = —alog(a) — (1 — a) log(l — a). (2.24) 


Finally the min entropy (n —» oo) is given by the inverse of the logarithm of the largest 
eigenvalue of pl 

S'|°°' ) = — log max [a, 1 — a]. (2.25) 

3 Generically the left-right decomposition of the operator is not as simple as in our example and depends 
on the structure of the Hilbert space of the CFT 
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Figure 1: 
in a state 


Finite value of the Renyi entropies as function of a for different replica numbers n 
excited by operator (2.14) 


All the Renyi entropies for n > 0 are sensitive to the special points a = 0 and a = 1 for which 
they vanish, reflecting the fact that for these values (2.21) becomes a product state. On the 
other hand, for a — f they saturate the inequality set by the Hartley entropy as expected 
from this maximally entangled mixture (see Fig [l] for comparisons). 


Let us stress that for general operators like EPR with different a each Renyi entropy gives 
a different constant dependent on a and only for a = 1/2 they coincide. This is the case 
studied in |9] and corresponds to the sigma operator in the Ising model. Hence, we expect 
that only for the rational CFTs all the Renyi entanglement entropies will increase by the 
same constant. 


Energy density 


The result (2.20) has been interpreted as the effect on entanglement due to an EPR pair 


whose members are receding from each other at the speed of light. To test this intuition, one 
could compute the time evolution for the expectation value of the energy density 

T tt (x,x) = -(T(x) +T(x)) , (2.26) 

in the excited state corresponding to O a |6j. The three-point function of the primary oper¬ 
ators with the stress tensor is universal and the final answer is given by 

,rr f u _ ( Ol(w 2 ,W 2 )Ttt(x,x)O a (w 1 ,w 1 )) 

\ tt\^^))O a , , _ W/-.f / - \\ 

{0 a {W l ,W 1 )0' a {W2,W 2 )) 

4 he 2 Ah e 2 

+ 


(Or + / - tf + e 2 ) 2 ' (Or + Z + f) 2 + e 2 ) 2 ' (2 ' 27) 

Indeed, the above expression describes two wave packets of width e propagating in the op¬ 
posite directions from the insertion point. The total energy injected to the system with the 
operator is equal to E = j7r(/ ' +/l) . 
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<r tt /h> 
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— t=0 

— t=3 

— t=8 


Figure 2: Evolution of the energy density for three different times. Plot for / = 0 and e = 1. 


The evolution of the energy density profile strongly supports the EPR interpretation of 
the operator excitation since the two wave packets resemble the pair (see Fig. [2]). Notice that 
the above computation is very robust and insensitive to the details of the primary operator, 


for example, the value of a in (2.14). 


3 Descendants examples 

In this section we compute the increase in the second Renyi entropy for few simple descen¬ 
dants. As we shall see, although some modifications on the energy density arise, the behavior 
of the Renyi entropy remains the same. 


3.1 Energy-momentum tensor 

Arguably, the most important non-primary operator in a CFT is the stress-energy tensor 
T(z). Recall that T(z) is a descendant of the identity operator with weight h = 2. It 
must be clear that the insertion of the identity operator into the system doesn’t cause any 
modification in the Renyi entropy. It is with the stress-energy tensor that we start our study 
of the effect of the insertion of descendants on entanglement. In the following, we compute 
the change of the second Renyi entropy Adue to the insertion of the chiral part of the 


stress-energy tensor T(z). For this case, equation (2.6) reads 


exp(-ASf^)) = 


(T(w 1 )T(w 2 )T(w 3 )T(w i ) } E2 


(3.1) 


«T( roi )T( W2 )) Ei r 

In order to compute the numerator of the above equation, we map the two-sheeted surface 
to the plane using the map (|2.7), under which T(z ) transforms as 


T( Wi ) = 


(1 


_2\2 r 


4 z?L 2 


T{Zi) + 


= ot{ [T(Zi) + 


(3.2) 












On the plane, the four-point function of the energy-momentum tensor is given by 


(T(zi)T(z 2 )T(z 3 )T(z 4 )) = (z 13 Z 24 ) 1 Gt{z) = (2132242(1 - z)) x Ft(z), 


(3.3) 


where fl8 


and 


g T (z) - j (1 + i + 


+ 2 c 


(!-*(!-*)) 
^ 2 ( 1-^) 2 ’ 


(3.4) 


T t {z) = - [z 4 (l - z) 4 + (1 - ^) 4 + 2 4 ] + 2c£ 2 (1 - ^) 2 [1 - z( 1 - z)} , (3.5) 

where the latter is constructed from Gt(z) by extracting the leading singular behavior. 


in terms of the cross ratios and read 


The contributions to (J3.l|) coming from the Jacobian and the denominator can be written 

TT 4 Qi- z 4 [\ — z') 4 

44>(3>’1 )TM>|, = c 2 /4 

Hence, we hnd that 


(3.6) 


exp (-A S^\z )) = 


^ 4 (1 - z) A 
c 2 /4 


Gt{z) + - + 


c 4 c 3 [1 — z( 1 — z)Y 


2 c 2 


16 4 z 2 ( 1 — z) 2 z( 1 — z) 


(3.7) 


where the terms in the bracket originate from the non-vanishing correlators of products of 

T{zi) + ft. 

Once again, as we take e —» 0, we either have z ~ 0 for t ^ [l, L + l] or z ~ 1 for 
t G [7,1/ + /]. Notice, that the pre-factor of (3.7) cancels the divergence of Gt(z) at those 
points. This is, moreover, the highest degree divergence inside the brackets. To make this 
fact manifest, we rewrite (3.7) as 


exp(-ASf(.-)) = ^ 


Ft(z) + z“( 1 - z)“ - + 


c 4 c 3 [1 — z{l — z)Y 


2 c 2 


Therefore, we hnd 


exp ( -A^ 2) ) ~ 


16 4 z 2 ( 1 — z) 2 z{l — z ) 


^t(O) t£[l,l + L] 

4 J- T ( 1) te[l,l + L] 


(3.8) 


From equation (3.5) we see that 


thus 


T{z —» 0) — T{z —>■ 1) — —, 


A .S';) 1 = 0 


(3.9) 

(3.10) 
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for all times. This way we confirm that both, primary I and descendant T(z) lead to the same 
behavior, A= 0, for any 1+ld CFT (not necessarily rational). As we will show, this 
phenomenon is not a coincidence and we will provide further evidence for it in the remainder 
of this work. 

A clarification at this point is in order. From the above formulas we can notice an order 
of limits conflict if we consider CFTs with large central charge. Namely in order to extract 
the interesting constant (in this case equal to zero) we first have to take e —> 0 and then 
c —» oo (for the opposite interesting order see App [E]) . 


3.2 Derivative of a primary 

In this section we consider a state locally excited by the first descendant of a primary operator 
O. We denote the descendant operator by 


dO(z, z ) = L_iO(z, z) = d z O(z, z). 


(3,11) 


Clearly, the conformal dimensions of dO are (h + 1, h), where ( h , h) are those corresponding 
to the primary. Let us start by computing the evolution of the energy density in states 
locally excited by (3.11). As mentioned above, the expectation value of the energy density 
is universal. As a matter of fact, we will only need the OPE between the stress tensor and 
the derivative of a primary field which reads [ 16] 


T(x)dO(zi,Zi 


2 hO(zi, Zi) (h + l)d Zi O(zi,Zi) &j.O(zi,Zi) 

X-Zi 


(3.12) 


(x - Zi ) 3 (x - Zi ) 2 

whereas the OPE with T(x) is the same as for the primary. 

The correlators with the stress tensor and descendants can be effectively written in terms 
a differential operator acting on the correlation functions of the primary operators (see also 
section 3.4). More precisely, we can write the three-point correlator as a sum over the residues 


from the OPEs of T and T with diOi, i — 1,2, taken inside the correlator and then pull the 
derivatives in front, such that the operator acts on the two-point function of the primaries 
O. Applying the differential operator and inserting the points (2.5), the expectation value 
of the energy density can be written as 


(T tt (x,x)) d0 = {T u (x,x))o + 


4e 2 [(1 -4 h)(x + l- tf + (1 + 4 h)e 2 } 


(2 h + 1) [(x + l - t) 2 + e 2 ] 3 
and the total energy injected into the system equals 

r f j it ( \\ 27 r(h + l + h) 

E ao = dx {Tu{x,x)} ao = -, 


(3.13) 


(3.14) 


as expected. Once again, the evolution of the energy density corresponds to that of two wave 
packets propagating in opposite directions from the insertion point. However, in the present 
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case there is an asymmetry in the amount of energy propagating on each side. Figure [3] 
clearly displays the asymmetry, where the right moving wave packet is seen to carry more 
energy. Naively one could expect that this different behavior of the energy density will imply 
the different constant increase in the Renyi entanglement entropies. We will show below that 
this turns out not to be correct. 


<7 tt /h> 


— t=0 

— t=3 

— t=8 


X 



Figure 3: Evolution of the energy density for three different times in the state excited by the 


descendant (3.11). For comparison, the dashed lines show the evolution in the state excited 


by a primary. Plot for l = 0 and e = 1. 


Let us first confront this with the left-right entanglement of the excited state. Formally, 


we can also think of the descendant operator (3.11) as acting on the left and right vacuum 
so that we have a new entangled stat^] 

d z O a |0) ® |0) = iV2a 

The normalized density matrix of the holomorphic or anti-holomorphic (left or right) movers 


a fa 


e zV2aA _ 


e -iV2a<t>\ 





(3 


for this state becomes (2.22) and therefore the corresponding Renyi entropies are again 


equal to (2.23). This suggests that, despite the different evolution of the energy density, 


the contribution to the entanglement Renyi entropy should be the same as for the primary 
operator. 

Now, we proceed to calculate the change in the second Renyi entropy wrt to the vacuum 
due to the insertion of the descendant operator dO. This can be obtained by computing 


exp 


-AS 


( 2 ), 


Z,Z = 


(dO(w 1 ,w 1 )dO(w 2 , w 2 )dO(w 3 , w 3 )dO(w 4 , w 4 )) s 2 
((dO(w 1 ,w 1 )dO(w 2 ,w 2 ))s 1 f 


(3.16) 


4 For simplicity we suppress the L/R indices that indicate left or right movers in the Minkowski signature 
or the functions of 2 or z in Euclidean. 
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The first step is to map the two-sheeted surface S 2 to the plane using the map (2.7). One 


must proceed with caution since descendant operators transform differently form primaries 
under conformal transformations (see App.B). For the first descendant of a primary we have 


dO(wi,Wi ) = (w() 




Wz 


dO(zi, Zi ) - h-jO(zi, Zi) 


w, 


= '' ,' 1 , [dO(zi, Zi) + pi 0{zii s;)] 


(3.17) 


Schematically, the correlator in the numerator of Eq. (3.16) can be written in terms of 
correlators on the plane as 


(dO(w l ,w 1 )...dO{w 4 ,w 4 )h 2 = ]Ja i a i ((dOdOdOdO) + ^pyOdOdOdO) + 

*=i ' j 


Y^PMoodOdO) + J2Wi(000d0) + A ( OOOO) 

j,k j,k,l 


i —1 


(3.18) 


The four-point function of the primary operator O is given by Eq. (2.16), whereas that 
for the first descendant takes the form@ 


(dO(z u z 1 )dO{z 2 , z 2 )dO(z 3 , z 3 )dO(z 4 , z A )) = (4 z A z 2 ) 


-2( ft + 1 ) (A~. = 


(4:ZiZ 2 ) ^Gdo(z,z). (3.19) 


The expression for Gao can be found in appendix [Cj equations (C.3) and (C.4). On the other 
hand, the two point function of 80 is simply 


{dO(wi,wi)dO(w 2 ,w 2 ))x 1 = - 2 h(l + 2h)w 1 2 l ' h+ 1 ' > w 1 2 h 
Now, we split the RHS of equation ( |3.16[ ) into two pieces 

exp z)^j = Mao (Gao(z, z) + E(z, z)) , 

where Mao is an overall pre-factor given by 

( 4 ^iz 2 )' 2(/l+1) (Aziz 2 y 2h nil OLiOLi 


Mao = 


,2 ’ 


((dO(w 1 , wi)dO(w 2 , w 2 ))sj 

and 5 contains all the terms of lower order in derivatives, i.e. 

E(z,z) = ( 4 Zl z 2 f {h+1) ( 4 z^f* (s^pyOdOdOdO) + YPjPk{°OdOdO) 

' j j,k 

4 

Y^PjPkPiiooodO) + n A (OOOO) 

j,k,l i= 1 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


5 Notice that this form assumes z 3 = —z\ and z A = — z 2 and for general points one cannot write the four 
point correlator of descendants this way. 
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To find the change of the second Renyi entropy, we use the Jacobian of the uniformization 
map (2.7) and the two point function (3.20) to show that the pre-factor can be written in 

(3.24) 


terms of the cross ratios as 


Mao = 


((1- z)z)^ 1+h \{l-z)zf h 


Therefore, we have 

exp (— AS^\z, z) 


Ah 2 (l + 2h) 2 
({l^z)z)^ 1+h \{l-z)z) 2 ~ h 


[' Gao(z,z ) + E(z,z)] 


4/r 2 (l + 2 h) 2 

For primary operators we have decomposed Go into 

Go(z, z) = (z( 1 - z))~ 2h (z( 1 - z))~ 2hjr o(z, z), 
where Tq is regular. Analogously, for the descendant we define 

Gao(z,z) = (z(l - z))~ 2[ - h+1 \z{l - z)y 2hjr ao(z,z), 


(3.25) 


(3.26) 


(3.27) 


where once again Tqo is regular and its relationship to Tq can be found in Eqs. (|C.7|) and 

Fao{z,z) + E(z,z) 


(C.8). Using these expressions, (3.25) becomes 

1 


exp(-ASf(2,z)) = 


4/r 2 (l + 2h) 2 


where we introduced 


l(z, z) = ((1 - z)z) 2{l+h \(l - z)z) 2 ~ h *(z, z) 


(3.28) 


(3.29) 


Now, we calculate A S^\z,z) after removing the cut-off. As discussed in section [ 2 ] we 
have to possibilities, either t E [/, L + 1} or t ^ [l,L + l\, for which the cross ratios become 
(z,z) — (1,0) and (z,z) — (0,0) respectively. Expanding 5 about the former yields 


<(z,z) ~(1 - z) 


Ah -2 


.(1 + 2 h) 2 
while for the latter we find 
Ah — 2 


*^£>( 1 , 0 ) + 


1 - 2 h 
1 + 2 h 


dfo( 1 , 0 ) 


+ zd? 0 { 1,0)+.... (3.30) 


l(z, z) ~z 


(1 + 2 h) 2 


^o(0,0) + 


1-2 h 
1 + 2 h 


dT 0 { 0 , 0 ) 


+ z dJ~o(0, 0) + 


(3.31) 


Thus, both of these contributions vanish. Finally, using Eq. (|C.7|) we find 

7(2) 


A SY ~ 


log T 0 (0,0) 
log ^o(l,0) 


t [/, + + /] 
t E [/, + + /] 


(3.32) 


which coincides with the result obtained for primary operators. The fact that we are consid¬ 
ering a descendant operator instead of a a primary has a tangible effect. Namely, the total 
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energy injected into the system is divided in different proportions amongst the left and right 
movers. In fact we would expect the action of each L_ n (L_ n ) to increase the height of the 
right (left) moving lump. Notice, however, that this unbalance doesn’t affect the jump in 
the second Renyi entropy of the system. It is only important that we have a pair of lumps 
propagating in the opposite directions and that one of the members of the pair is inside the 
entangling region of an arbitrary shape (in Id either a finite or semi-infinite interval) while 
the other member remains outside. This strongly hints to the topological nature of this 
quantity (see also [25] and the discussion section). 


3.3 Both derivatives 

In this section, we briefly consider the change in the second Renyi entropy due to the insertion 
of the second order descendant 


ddO(z, z ) = L_iL_iO(zi z) = d z d z O(z, z). 


(3.33) 


The energy density profile corresponding to this descendant is given by 

IT ( \\- —It( \\ , ^e 2 [(1 — 4h)(x + / — t) 2 + (1 + 4h)e 2 ] 

\Ttt{x, x)) aao —(T tt {x,x))o H- 


+ 


(2 h + l)[(x + 1 — t) 2 + e 2 ] 3 
4e 2 [(1 — 4 h)(x + l + t) 2 + (1 + 4 h)e 2 ] 


(2 h + 1) [(x + l + t) 2 + e 2 ] 


21 3 


(3.34) 


and its time evolution is depicted in Fig. [4] Moreover, the total energy injected into the 
system is 

2n(h + h + 2) 


EddO ~ 


(3.35) 


Notice that the left/right symmetry of the energy carried by the wave packets is restored. 

The task of calculating the A S y A corresponding to the insertion of ddO is analogous to 
those performed in the previous sections. Nevertheless, there are some complications that 
we wish to point out. First, the transformation of the operator ddO under conformal maps 
is more complicated, it contains contributions proportional to ddO , dO and O. The explicit 
form of this transformation can be found in Eq. (B.2). Once again, we find an expression of 
the form 


exp (-A^ 2) (z, z)) = A r Beo [ Q~ dao {z , z) + ... ] , (3.36) 

where the dots contain the elements of lower order in derivatives. Furthermore, the normal¬ 
ization factor reads 


/((l - ^) 2 ( 1+/l )\ (((1 - z)z) 2 ^ 1+h )\ 

V 4h 2 (l + 2h) 2 ) ^ 4h 2 (l + 2h) 2 ) ' 


(3.37) 
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<T tt /h> 



Figure 4: Evolution of the energy density for three different times in the state excited by 
ddO. The dashed lines show the evolution in the state excited by a primary. Plot for l = 0 
and e = 1. 


As well, the function Qqqq becomes more convoluted; nevertheless, it can be expressed in 
terms Q 0 , see Eq. (C.10). Despite these complications, we can show that 

exp (—AS^\z,z)j — (4:h(l+ 2h)h(l + 2h)) 2 [ Fqq 0 {z,z), (3.38) 


where Tqqq is given by (C.12). Finally, as we remove the UV cut-off we End that 


A Sf ~ 


-log Jo (0,0) 
-log Jo (1,0) 


t <f [l, L + l] 
t £ [/, L T /] 


(3.39) 


once more. This clearly shows that the propagation of the energy in the state excited by the 
descendants only provides a qualitative support for the quasi-particle picture. 


3.4 Descendants at level 2 

It is always possible to write the correlators of descendants in terms of the action of certain 
differential operators on the correlator of the primary operators (family’s parent). Whenever 
the correlator in question involves only (powers of) L_i and L_i the Virasoro generators can 
be pulled out of the correlator as simple derivatives. However, if there is any descendant 
constructed with an L_ k with k > 1 the situation changes. For example, if there is one such 
descendant in the correlator then P3] 

(0(wi,wi).. .L_ k O(wi,Wi)... 0(w N ,w N )) = £ M k (0(w 1 ,ui 1 )... 0(w N ,w N )), (3.40) 


where 


Ai) _ 

J ~'-k ~ 



h{k — 1 ) 

k 

w■■ 

31 



(3.41) 
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However, if the correlator contains more than one descendant the appropriate commutation 
of the Virasoro algebras at different points must be taken into account. This is manifest 
already for level k — 2 and we discuss it below. 

(o') 

In the following we compute A S y A due to the insertion of the operator 

0 ( ~ 2 \z,z) = L_ 2 0(z,z). (3.42) 


The crucial correlation function that will contain the relevant constant is the four point 
function of 0^~ 2 \ which appears in the numerator of Eq. (2.61 6 We proceed in a standard 
way [16], first we write the correlator as 


(L^OiO^O^O^ 


> = £ 


i =2 J C(m) 
and then make use of the OPE 


d z 

—X Zl - z )-'(0 1 ...(T(z)0, 


(- 2 ) 


(3.43) 


T(z)Qt 


-2), 


W) 


(f+ 4 h) 0(w) + 3 dO{w) + (h + 2)0 ( ~ 2 \w) + dO^- 2 \w) 


(z — wy 


(z — wy 


(z — w) r - 


w 


(3.44) 


Inserting this OPE into (3.43) and calculating the residues, we can express the four point 
function in terms of correlators containing one primary and three descendants. Repeating 
this procedure, it is possible to reduce the number of descendants in the RHS until the full 
correlator is written in terms of the four-point function of the primary operators. The full 
procedure is rather cumbersome and the result is given by (D.6). 

Once more, after insertion z 3 = —z\ and z 4 = — z 2 , we define 


with O t = 0(zi, Zi ), as well as 


exp ( 

where the pre-factor reads 


(42i2 2 ) 2(h+l) (4Z!Z 2 ) 2h C/ o( -2)(z,z), 

(3.45) 

r 2(h+1) ( z(i-z))~ 2h j- 0( - 2) (z,z ). 

(3.46) 

y can be written as 


—A/"o(- 2 ) [ Qq(- 2 ) (z,z) + ...] , 

(3.47) 


(4zi2 2 ) 


-2(h+2) 


{^ZiZ 2 ) 2I> n?=l a i a i 


° } {{0(- 2 \w 1 ,w 1 )0(- 2 \w 2 ,w 2 ))x i y 

Using the OPE ( |3.44 ), we can show that the two-point function is given by 

{0^ 2 \w 1 ,w 1 )0 { ~ 2 \w 2 ,w 2 )) Tj1 = ^ (c + 2h(9h + 22))w^ 2 2(h+2) w^ 2 2h . 


(3.48) 


(3.49) 


6 In the limit of (z, z) —> (1,0) all the other correlators are suppressed by powers of (1 — z) 
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Finally, we extract T 0 (- 2 ) from Eq. (D.6), and as we remove the UV cut-off e find 


-p ( -\ 1 , . ohfnh oo'o 2 / ^ 0) t £ [l,L + l] 

^o(- 2 )(z,z) ~ 4 (c + 2h(9/i + 22)) | jr o ( 1)0 ) ie[Z,L + Z] 


(3.50) 


Therefore, we find that for the insertion of the operator Ch 2 ) we have 


a r.(2) ^ f ~ ^0(0, 0) t£[l,L + l\ 

! 1 — log J r o( 1,0) t £ [/, T + /] 


(3.51) 


as expected. 

Clearly the procedure becomes more involved once the level of the descendant increases 
but the increase in the Renyi entropy remains fixed by the two limits of the original J r o{ z , z ) 
from the correlator of the primary operators. 


4 General derivation 


Based on the explicit examples in the previous sections, we can clearly deduce that the 
relevant contribution to the Renyi entanglement entropies is the same for the whole conformal 
family. We show it for the second Renyi entropy but the generalization to higher n should 
be straightforward (though very involved computationally). 

The technical reason for the constant is the following^ On one hand, under conformal 
transformation, every descendant of level (n, ft) transforms as 

z) = (f) h+n (/')“ [©<"'*>(/, /) + ...] (4.1) 


where the ellipsis stand operators of the lower dimension that in the correlation functions 
give rise to lower order singularity. This way, once we employ the map (2.7), the Jacobian 
from (4.1) with the square of the two-point function on £1 (2.6) combine into 


P = - z)f +n [z(l - z)\ 


—\lh-\-n 


(4.2) 


with M 2 being the constant in the two-point function of the descendants. 

When we extract the constant contribution to the entanglement Renyi entropy we take 
the limit of (z,z) —> (1,0). The only terms that survive this regime must cancel the above 
pre-factor and are confined to the correlation function of (n, n) descendants only. 

On the other hand, the correlators of descendants on the plane can be obtained from the 
correlators of the primary operators that we define as [^] 

(O a (zi,zi)Ol t (z2,Z2)O a (z 3 ,z 3 )Ol t (z4,Z4)) = \z 13 z 2 4\~ Ah \z(l - z) |~ 4 / ' T a {z, z), (4.3) 

7 Similar logic has been applied in )24j 

8 This is a form that we assume with F(z, z) being smooth at 2 —» 1 as well as z —> 0. This is holds for all 
the cases that we are aware of. 
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by acting with a (complicated but in principle straightforward to derive) differential operator. 
In the limit of (z,z) —» (1, 0) the only contribution that has enough singularity to cancel the 
above-mentioned pre-factor comes from the terms in which the differential operator leaves 
T a (z,z) intact and contains precisely the square of the norm Af 0 ( n ,n) of the two-point func- 
tioiiB 

This way the non-trivial constant obtained from the correlators in the limit of (z, z) —> (1, 0) 
is the same for all the members of a given conformal family. By the same token, one can 
show that as ( z , z) (0, 0) the Renyi entanglement entropy is unchanged. 


5 Conclusions and discussion 

In this work we have demonstrated that, in 1+1 dimensional CFTs with Virasoro symme¬ 
try, local excitations by operators in the same conformal family increase the second Renyi 
entanglement entropy by the same constant. For rational conformal held theories this entan¬ 
glement constant is the logarithm of the quantum dimension (obtained in |9j) of the primary 
operator that represents the family. We have checked this on explicit examples and outlined 
a general derivation. A generalization to higher Renyi entropies should be also possible and 
we expect it to work in exactly the same way as for the second Renyi entropy^ 

Our results strongly support the existence of universal features of entanglement for locally 
excited states. So far only the entanglement entropy of a block in the ground state of a local 
hamiltonian has been used to fix the central charge of the CFT that govern the critical points. 
It would be interesting to perform numerical study of the evolution of the entanglement 
entropy in a chain (e.g. Ising chain) excited by a local operator. Similarly to [32J, by tuning 
the parameters close to criticality^ it should be possible to verify the CFT prediction for 
the constant contribution to the entropy [33] . 

Let us also discuss two different perspectives that we hope can lead to a simpler proof of 
our result and shed more light on the physical meaning of the constant contribution to the 
Renyi entanglement entropies. 

We have stressed that the constant contribution to the Renyi entropies is equivalent to 
the entropy between Left and Right (chiral and anti-chiral) movers in a quantum state [3] 

\i, LR ) = O(z,z)\0) L \0) R (5.1) 

9 We have verified it up to level 2 and for arbitrary powers of derivatives 

10 Even though general 2n-point correlators of descendants will have a very complicated structure, the late 
time behaviour of the invariant cross-ratios is universal. Therefore, at least in RCFTs, one could possibly 
use the factorisation and the fusion transformation for any members of a conformal family to formally 
demonstrate our claims for general n. The same should be true for n < 1 and it would be interesting to 
check it for the logarithmic negativity (see e.g. mm 
n Generically, local operators on the chain, at the critical point, will correspond to primaries, descendants 
or the linear combinations of thereof. 
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Formally, this means that we can write the state in the Schmidt form 


V!>lr) = \ir) 

i 


with Schmidt bases | %l) for L and | in) R and Schmidt coefficients satisfying 




(5.2) 


(5.3) 


where number of non-zero piS is called the Schmidt number of state | i/ilr) [15] - This way the 
constant contribution to the n-th Renyi entropy can be computed as 


A S {n) 



(5.4) 


and comparing (5.4) and (2.12) we get for rational CFTs: ^oo[o] = Ei Pi 


The fact that this number is the same for all the members of the family of a give descen¬ 
dant indicates that only the elements of the Schmidt basis | ii) (|*r)) can change from one 
descendant to another however the Schmidt coefficients remain the same. 


In principle one should be able to express the excited state (5.1) in the left-right or¬ 
thonormal basis of descendants and then transform to the Schmidt basis using the singular 
value decomposition. This can be done for some excited states that have a natural left-right 
decomposition (see [201 ED [221 [23]) but it would be interesting to preform it for ( |5.1[ ) even 
in known RCFTs. 


Let us also point that the logarithm of quantum dimension appears naturally in the 
computation of topological entanglement entropy in 2+1 dimensions. In fact it was shown 
in ra that for a given region A the increase of the topological entanglement entropy due to 
excitation a is equal 

A S top = log (da) (5.5) 

where d a is the quantum dimension of the excitation. 


Even though rational CFT are naturally linked to topological field theories in 2+1 di¬ 
mensions (e.g. FQHe), the precise link to the 2d CFT technology employed above is not 
obvious. However, based on the connection between the topological entropy and the bound¬ 
ary entropy [2Ej, or left-right entropy and TQFT in 2+1 [22J 12 it is not unlikely that the 


constant contribution to Renyi entanglement entropies in 1+1 d CFT is equivalent to the 
change in topological entanglement entropy in 2+ld. 

If this happens to be the case, our results for descendant^] imply that the contribution from 
the whole topological sector (conformal family) is the same. It would be very interesting to 
sharpen this connection and we leave it for the future work. 


12 See also [28] for an interesting application to black holes in AdS^ 
13 That give rise to wave functions for excited states (see e.g.^) 
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A Correlators of the EPR operators 


We have used the EPR operators 

O a (zi, Zi) = V^e iV * a * {zi '* i] ± = y^Oj +) ± VT^O G) 

O f a ( Zi , zt) = ± Vl^ae iV2a ^ Zi ^ = Va Ol~ ] ± Vl^i,0 G) (A.l) 

with conformal dimension h = h = a 2 and parameter a G [0,1]. 

Here we compute the four-point correlator on the complex plane 

G o a = {O a (zi, zi)Ol(z 2 , z 2 )O a (z 3 , z 3 )Ol(z 4 , z A )) (A.2) 

Since the non zero correlators of the vertex operators are the ones that have the zero sum of 
the operators dimensions 

(JJ e iV2aG{Zi ’~ Zi) ) =Y[\ Zi - Zj | 4a<a U J2 ai = 0 A- 3 ) 

i i<j i 

we denote 

^(ai,02,03,04) = (oS ai) of 2) 0^ a3) Oi a4) ), a.i e {+, -} (A.4) 

And the non-zero contributions to the correlator comes from 

Go a = (a 2 + (1 - a) 2 )G (+ ’ _,+,_) + 2a(l - a) (G (+,+, “ + G (+ ’“ ~ +) ) (A.5) 

where we used 


G (+, = G (-,+,-,+) j q (+,+, 

The remaining correlators are 


-) = g g '-' + ’ + \ = G ( -’ + ’ + ’- } (A.6) 


|U3^24| 4h |^(l - z)\ Ah 


.-) = . 


\z\ 8h 


\zi 3 Z 2 4 

4 h\ 

\z(l — z)\ Ah 

,+) = . 

| 

1 - 

00 

Gt 

1 


^13^24 

4 hi 

\z(l — z)\ Ah 


(A.7) 
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The four-point correlator of such operators on the complex plane is given by 


G o a ~ 1^13^241 U ‘G(z,z) (A.8) 

with the canonical 4-pt function 

G{z,z) = (O a \O a (z,z)O a (l,l)\O a ) = \z(l-z)\~ Ah R{z,z) (A.9) 

where 

7{z, z) = a 2 + { 1 - af + 2a(l - a) (|z| 8/l + |1 - z| 8h ) . (A.10) 


B Transformation of descendants 


The transformation of descendant fields (of a primary </> with dimension h ) under finite 
conformal map f(z) has been worked out in [T3] (see also [S]). For the first few descendants 
it can be obtained from 

G ( ~ k \z) = e RoLo (l + R.L, + -R\L\ + R 2 L 2 + 0 ( " fc) (/) 

where 

pi p 1 w f" p Sf 

Ro = log / , Ri = -Rq = A 2 = — 

where Sf is the Schwarzian derivative of f. Using the Virasoro algebra and the commutation 
relations of L’s with the operator modes (f) n we can derive 


<t>(z) 

4> { ~ l \z) 

^~ l \z) 

G ( ~ 2 \z) 




(f ) h+1 

/\h-\- 2 




(/ 


+ (2 h + 1 )jr 2 G { - l \f) + h ({2h + 1)^ + p j m 


(/') 


f\h-\- 2 


^ 2) (/) + 1} (/) + (^r + (4/* + 


4/" 


57 

6/ /2 


0 (/) 


(B.l) 


The first two descendants correspond to derivatives and the transformation rule can be also 
reproduced from 


= o, !(/')U(/)l , =« f/oUa) 


It is also easy to check that for h — 0 (the identity) the transformation of 2 ) (z) reproduces 
that standard transformation for T(z). 
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We have also used the transformations 


dd<p(z, z) = (f)“ +1 (/') s+1 sa<Kf, D + D + /) + hh^ 2 f) 


dd<Hf,f) = (f‘) {h+1) (f') ( ' 1+1) dd<t>(z,z)~hyd</>(z,z)-hL-d<t>(z,z) + hhyL-<l>(z,z) 


(B.2) 


C Correlators at level 1 


In the following, we display some results concerning the four-point function of the operator 
dO on the plane. If we define 

(0(zi, zi)0(z 2 , z 2 )0(z 3 , z 3 )0(z 4 , z 4 )) = (z 13 z 24 ) 2h (-13-24) Go( z > z), (C.l) 


as well as 14 


{dO(zi, Zi)dO(z 2 , z 2 )dO(z 3 , z 3 )dO(z 4 , z 4 )) = (-13-24) 2(h+1) (-13-24) 2h Gao(z, z ), (C.2) 


we can show that Go(z,z) and Gao{z, z) are related by 


4 


Gao(z,z ) = ^9i{z) d l Go(z,z ) 

1=0 

(C.3) 

where 


9o{z) = 4/i 2 (l + 2 h ) 2 , 

9 i(z) — 2(1 + 2h)(2z — 1) , 
g 2 (z) = 2 + 2 z(Ah 2 -7- 2/t)(l - z ). 
g 3 (z) = —4 2(1 - z)(2z - 1) , 
g 4 {z) = (1 - -) 2 - 2 ■ 

(C.4) 

Moreover, if we define the functions 


To{z, z) = (z( 1 - z)) 2ft (S(l - 2)) 2S 5o(z, z) , 

(C.5) 

and 

T d0 {z,z) = (z( 1 - -)) 2(/l+1) (5(l - z)) 2h Gao(z,z) , 

(C.6) 

we can show that 

4 


Fdo{z,z) = ^2fi(z) d l Td(z, z) , 

1=0 

(C.7) 

14 This form assumes Z 3 = —24 and 24 = — 22 
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where 


fo(z) =4h 2 [(l + (-1 + z)z) 2 + 4/? 2 (1 + 3(—1 + z)z ) 2 
—4/i(—1 + (-1 + z)z (-2 + 5(—1 + z)z))\ 
fi(z) = — 2(—1 + z)z(--l + 2z)(z + 10h(—l + z)z — z 2 
— 40h 2 (—1 + 2)2 + 16/r 3 (l + 3(—1 + z)z)) 
f 2 ^z ) = 2 ( —1 + ^) 2 z 2 (l + 7(—1 + z)z + /i(—6 — 34 (—1 + z)z ) 

+ 4h 2 (3 + 11 (—1 + z)z)) 
f 3 (z ) = - 4(—1 + 2/i) (-1 + zfz 3 {-l + 22 ) 

/ 4 (z)=(-l + z)V. (C.8) 


Note: At the points (z, z ) —> (0,0), (1, 0), the only contribution to Tqo{z , 0) that survives 


comes from the fo(z) term in (C.7), which in terms of the coefficients of (C.3) reads 


Fdo{ z i z ) 2h(l + 2 h) [g 2 ( z ) + 2(1 + h)(( 3 + 2/1)74(2) — 73(0))] J 7 ’0(0, 0), 

where £3(0) = 0(1 - 0)73(2) and g 4 (z) = 0 2 (1 - z) 2 ^ 4 (z). 

On the other hand, for the descendant BdO we have 


(C.9) 


where 


Gdoo{ z , z ) = 9 kA z ^ z ) d k d l g 0 ( z , z ) 

k,l =0 

9k,i{ z i z ) = 9k{z)gi{z), 


where the functions gi(z) are defined in (C.4). 


(C.10) 


(C.ll) 


? 8 do(z, z ) = {A 1 - 0)) 2(fe+1) (-( 1 - z )) 2 {h+ 1 ] Gddoi z i z ) , 


(C.12) 


D Correlators at level 2 


In this appendix we present some details of the computation of the correlation functions of 
descendants 0^ 2 \w,w). The key ingredient (for any operator) is the OPE with the energy 
momentum tensor that in this case reads 


~ + , 39QM (ft + 2)0<- 2 >M SO^Hw) 

(0 — w ) 4 (z — w ) 3 (z — w) 2 z — w 

Denoting 0^ 2 \wi) = 0-~ and 0(wi ) = 0 % we compute on the complex plane 


(D.l) 


(0[~ 2) 0 { 2 ~ 2) ) = 


£ + 4 h 3d 2 

+ -T + 

™12 


™t 2 


(h + 2) d 2 

2 ' 

™12 


W 12 


h d x 
W\2 W 12 


(Oi0 2 ) 


= - (c + 2h(9h + 22)) w 12 


- 2 (h+ 2 ) 


(D.2) 
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where Wij = Wi — Wj. The products of operators should be understood as acting from left 
to right on the correlators. When taking the contribution from the residues (see e.g. 6.6 in 
[S]) there are always few ways to start evaluating the correlator and we use the convention 
from left to right as well. 


Analogously the three-point function can be written 

(i o[- 2) o { ~ 2) oi- 2) ) = [p 1>2 £® + Pi, 3+ ^J 1 1 ) ) p 2 ,3 + (o^o ,), 
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Finally, the four point function can be written as 

(0[~ 2 ) 0 { 2 - 2 ) 0t 2 ) 0[- 2) ) =\v li2 fv 3 , 4 + lff 2) £ty +p 1)3 fv 2A + lff 3) £W 


+Pl,4 ( P 2 ,3 + Pp,4)^- 2 ) + ^(ij ( ^A,2^-2 + ^4,3^ 
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•(3) 


,( 2 ) 
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i -iy(l)7"(4) 
Xn (l) X (l,4) 


-r) I 7 -C 2 ) <^(3) 

2V 2 ,3 + X (l,4) X '-2 




where 


~ ~7.) ^ hk = hi = h and h p = h + 2 for p ^ {fc, /}. 






(D.3) 

(D.4) 

(D.5) 


(D.6) 


(D.7) 


One can verify that setting h = 0 in the above operator reproduces the four-point function 

(pk. 


E Comments on large c and holography 


The late time value of the Renyi entanglement entropies is non-perturbative. To illustrate 
this fact we consider again the state excited by the descendant of identity T(z). In the 
computation of the late time value (semi-inihnite interval) of the second Renyi entanglement 
entropy we use the (3.7). If we first take the limit of large central charge, the contribution 
to the second Renyi entropy becomes 


zt c 

AS'( 2 ' ) ~81og-2 log- (E.l) 

6 Z 

This is the logarithmic growth with time obtained in [6j (the negative part comes form the 
normalisation) and it should be possible to perform a similar computation for the entangle¬ 
ment entropy and match it with gravity using |33J as in (7J. 

On the other hand, as we saw in the text, if we know the answer for the full correlator 
(all quantum gravity effects taken into account) we can first take e —)• 0 and obtain the 
entanglement constant. 
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